Thermodynamic formalism for the Lorentz gas with open boundaries in d dimensions 
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A Lorentz gas may be defined as a system of fixed dispersing scatterers, with a single light particle moving 
among these and making specular collisions on encounters with the scatterers. For a dilute Lorentz gas with 
open boundaries in d dimensions we relate the thermodynamic formalism to a random flight problem. Using 
this representation we analytically calculate the central quantity within this formalism, the topological pressure, 
as a function of system size and a temperature-like parameter j3. The topological pressure is given as the sum 
of the topological pressure for the closed system and a diffusion term with a /9-dependent diffusion coefficient. 
From the topological pressure we obtain the Kolmogorov-Sinai entropy on the repeller, the topological entropy, 
and the partial information dimension. 

PACS numbers: 05.45.-a, 05.70.Ln, 05.90.+m 



INTRODUCTION 

In the seventies of the last century Sinai, Ruelle, and Bowen 
developed a formalism for calculating dynamical properties 
of chaotic dynamical systems UL^y]. This formalism re- 
sembles very much Gibbs ensemble theory and therefore was 
given the name thermodynamic formalism. From a partition 
function defined as an integral over phase space of a particular 
weight function, one may derive several dynamical quantities, 
such as the Kolmogorov-Sinai (KS) entropy, the topological 
entropy, the escape rate, or the partial information dimension 

Although the power of this formalism has been widely rec- 
ognized, there are only few examples so far of physically in- 
teresting systems for which the topological pressure and re- 
lated functions could be evaluated explicitly. One such sys- 
tem is the dilute random Lorentz gas. Generally spoken a 
Lorentz gas is a system consisting of fixed dispersing scat- 
terers, among which a single light particle (or alternatively 
a cloud of mutually non-interacting light particles) moves 
about, making specular collisions on each encounter with a 
scatterer. The scatterers may either be placed on a periodic 
lattice or at random positions in space. Usually the scatterers 
are taken to be disks, in d — 2, spheres, in d = 3, or hyper- 
spheres, for d > 3. If these are placed on a periodic lattice the 
resulting system is also known as the Sinai billiard |6i]. In fact 
this is also the model Lorentz 1 13] had in mind originally for 
describing the transport of conduction electrons in metals (still 
in the pre-quantum era). However the kinetic equation he pro- 
posed to describe this system, presently known as the Lorentz - 
Boltzmann equation (see for instance OlllBlX m f act is more 
appropriate for the model with random scatterer locations, at 
low density of scatterers. The Lorentz-Boltzmann equation 
and some of its generalizations to higher densities I KIM KM 211 
allow for analytic calculations of transport coefficients and 
other fundamental nonequilibrium properties. 

Recently, several dynamical properties have been calcu- 
lated analytically for open and closed random Lorentz gases 
by using an extended Lorentz-Boltzmann equation approach 



(HUl Van Bei J eren and Dorfman (H alternatively 
used the thermodynamic formalism to calculate the topologi- 
cal pressure and related quantities for a dilute d-dimensional 
random Lorentz gas in equilibrium. For 2 dimensions this was 
extended by the present authors IT^l to the case of a dilute ran- 
dom Lorentz gas under the combined actions of a driving field 
and an isokinetic gaussian thermostat. 

THERMODYNAMIC FORMALISM 

In chaos theory a central role is played by the time evolu- 
tion of infinitesimal volumes in phase space. For chaotic sys- 
tems such a volume element will grow in some directions and 
shrink in others. The factor by which the projection of the vol- 
ume element onto its unstable (expanding) manifold increases 
over a time t is called the stretching factor. For an infinitesi- 
mal volume centered originally around the phase point (r,v) 
we will denote its value as A(r, v, t). Similarly, the contrac- 
tion factor is the factor by which the projection onto the stable 
(contracting) manifold decreases over time t. In analogy to the 
Gibbs ensemble formalism of statistical mechanics, we may 
define a dynamical partition function Z(J3,t) as a weighted 
integral over phase space. For closed systems, indicated by a 
subscript 0, points in phase space are weighted by powers of 
their stretching factor, according to 



dfi(r, v) [A(r, v, t) 



1-/3 



(1) 



where the integration is over an appropriate stationary mea- 
sure. In systems with escape, phase space trajectories are re- 
moved from the ensemble if they hit an absorbing boundary. 
In this case the definition of the dynamical partition function 
has to be generalized to 

Z(J3, t) = J dfi{r, v) [A(r, v, t)} 1 ^ 3 XtP (r, v, t), (2) 

with xtr{r,v : f) — if the trajectory starting from (r,v) 
at time hits the absorbing boundary before time t and 
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Xtr(i",v,t) = 1 otherwise. In our analogy to statistical me- 
chanics the parameter (3 behaves like an inverse temperature 
and as the analogon of the Helmholtz free energy we obtain 
the topological pressure P{0) as 



P(f3)= lim jln Z(p,t). 



(3) 



The dynamical entropy function h((3) is given by the Leg- 
endre transform of P((3) as, 



h(0) = P(/3) - /3P'(/3), 



(4) 



where P'(f3) — dP((3)/d(3. For special values of (3 
the dynamical entropy function can be identified with dy- 
namical properties, because for long times we have A ~ 
exp(i J^t ^')> where J^t ^ * s tne sum of all positive Lya- 
punov exponents A^. Specifically, h((3) equals the topological 
entropy h top for (3 = and the KS entropy h^s f° r /3 = 1. If 
the system is closed, the KS entropy equals the sum of positive 
Lyapunov exponents because Pq(1) = 0, as follows directly 
from Eqs. and (|3}. 

However, for open systems P(l) = —7, where 7 is the 
asymptotic escape rate. The relationship P'(l) = \ 
still holds, but now the Lyapunov exponents are defined on 
the repeller, i.e. the subset of phase space from which no tra- 
jectories escape. The point where P((3) intersects the /3-axis 
can be related to the partial Hausdorff dimension, which is a 
fractal dimension of a line across the stable manifold of the 



attractor. Another fractal dimension associated with the topo- 
logical pressure is the partial information dimension which is 
given by the intersection point with the /3-axis of the tangent 
to P(/3) at P(l). For the closed system both dimensions co- 
incide and are equal to 1. 



DYNAMICAL PARTITION FUNCTION FOR THE OPEN 
LORENTZ GAS IN d DIMENSIONS 



In this section we will calculate Z((3,t) for a dilute ran- 
dom open Lorentz gas in d dimensions, with N fixed (hy- 
per)spherical scatterers of radius a distributed randomly in- 
side a finite volume V. In addition there is one point particle 
moving with velocity v along straight lines between specu- 
lar collisions with the scatterers. If the particle leaves V it 
escapes. "Dilute" implies the condition na d <C 1, with the 
density n = N/V. 

The condition of diluteness allows us to approximate the 
dynamics of the light particle as a random flight, in which each 
trajectory between subsequent collisions is sampled from an 
exponential distribution of free path lengths and the collision 
parameters of each collision are sampled from a distribution 
corresponding to a homogeneous bundle hitting the scatterer. 
In other words: all effects resulting from recollisions are com- 
pletely ignored. Under these approximations the dynamical 
partition function may be rewritten as 



dr / dvS(\v\ - v ) dtf-dtt dp,{a x ) ■ ■ ■ dp(m) 

J l=0 Jo J 



6 ( t - ^ ti exp 



i=l 



K + Po(/3))(i-5> 



f 1 \ 1 

n + v[{t - ]T U) [] W{a ,t 3 ) X (r ). (5) 

v z=l / j=l 



Here Pq is the topological pressure for the closed Lorentz gas 
at the same density in equilibrium; a. L denotes the collision 
normal at the ith collision and d^,(a) denotes the probability 
measure for scattering with collision normal a; <d(x) denotes 
the unit step function, i.e. Q(x) = 1 for x > and Q(x) =0 
for x < 0; Vd = 1 / fd is the average collision rate, given for 
dimension d as 1 141 
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d-l„ 



(d-i)r(^i) 



(6) 



Tj is the position of the light particle at the jth collision, 
and x( r ) i s m e characteristic function satisfying x( r ) = 1 
for r inside V and otherwise. We implicitly assumed that 
the boundary is not concave, so if both x( r j) an d x( r j+i) 
equal unity the same is true for the characteristic function of 
all points in between. In addition, v[ is the velocity of the 



light particle just after the Zth collision. Finally, W(a, t) is an 
effective free flight transition rate, defined as 

W(a, t) = Vde -&W+ V *'> t [A d (9, t)] 1 ' . (7) 

Here the stretching factor is given by fill [nil 

A d (6 i t)=(^\ \cos9) d ~\ (8) 

where 9 is the scattering angle, defined through cos 9 — —v-a, 
with v the unit vector along the velocity of the light particle 
before the collision, see Fig.^ Since in Eq. the only de- 
pendence of the integrand on dj occurs through cos 6^ , the 
integrations over dp(dj) may be reduced as 

, ,7T/2 

J dn{dj) -> (d - 1) J d9j cos 9 3 sin d " 2 9 r 
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FIG. 1 : Sketch of a trajectory in two dimensions of the light particle 
starting form the initial point and following the first two collisions. 



Note that in Eq. the actual free flight distribution has been 
changed to an effective free flight distribution by multiplica- 
tion by the (1 — /3)th power of the stretching factor and by the 
factor exp(— Po(/3)t). Similarly, the distribution of collision 
normals has been changed to an effective distribution. Indeed, 
after this rearrangement the integral of W(a, t) over dfi(a) 
and positive t equals unity 1 14]. The second moments of both 
time and displacement for the effective distribution W(a, t) 
are well-defined, and the resulting effective random flight pro- 
cess for given initial direction gives rise to a convergent aver- 
age displacement after n steps in the limit n — > oo. Therefore, 
on large time and length scales it leads to a normal diffusion 
process, with a /3-dependent diffusion coefficient D((3). 

On division of the logarithm of Eq. (0 by t the first factor 
just simply reduces to Po((3). The contribution from the re- 
maining factor may be interpreted, up to a normalization fac- 
tor, as the average survival probability of an initially homoge- 
neously distributed cloud of light particles of equal energy, un- 
der the effective random flight process with absorbing bound- 
aries described above. Since on large time and length scales 
this becomes an effective diffusion process, this survival prob- 
ability will behave as exp(— fcgZ?(/3)i) for long times, with fc 
determined by the geometry of the system and to some ex- 
tent by the boundary conditions on the random flight process. 
Therefore the second contribution is of the form — fegD(/3). 
Combination of these two brings us to the main result of this 
paper: the topological pressure for the dilute open random 
Lorentz gas may be expressed as 



P(/3) = P ((3) 



k 2 D(/3). 



(9) 



What remains to be done now is finding an explicit expres- 
sion for D((3) as a function of both (3 and d. This task is not 
particularly hard. The Green-Kubo expression for the diffu- 
sion coefficient as a time integral of the velocity autocorrela- 



tion function gives rise to 
1 



dt (v(0) ■ v(t)) 



v 
~d 



TO 



2 = 1 



[cos(ir - 20)) 1 



(10) 



The first term is the contribution from free flights from the 
initial time until the first collision, the next terms result from 
free flights between the Ith and (7 + l)th collisions. Since 
the collision cross-sections are isotropic and all collisions are 
uncorrected, the average direction of the velocity after the Zth 
collision is proportional to (cos(7T — 20))'. Then we obtain 
from Eqs. (|5j and @ the time averages 

/n-2\ 



1 



TO 



2(r 



(r) 



1 



+■ W) 

together with the angular average 

(cos(tt - 26>)) 



1 



(d -!)(!-/?) 



-0(3 - d) 



(ID 
(12) 

(13) 

2(d-l) + /3(3-d) 

To understand the first equality in Eq. ( II 11 one should real- 
ize that the probability of finding the initial time on a free- 
time stretch of length r is proportional to r and that the av- 
erage time until the first collision under these conditions is 
just t 1 2. Furthermore Eqs. (II It and ill 2b only make sense for 
(3 < d/(d—l) since for larger /3-values the effective free-time 
distribution shows a divergence, see Eqs. and (|8j. Simi- 
larly Eq. d 1 3I> has to be restricted to the range (3 > — 1, for 
d = 2, respectively (3 < (d - l)/(d - 3) for d > 4 (for d = 3 
there are no restrictions). 

The topological pressure for the closed Lorentz gas in equi- 
librium is obtained from the first pole of the Laplace transform 
of the dynamical partition function il4lfl9il and for the closed 
Lorentz gas in equilibrium it reads 

(d-l)(l-/3) 

T (d + j3 - d/3) 



2v 



;^)r(— 



-l+/3(3-ri) 
2 



r ( d - 1 + 



>(14) 



Note that for the closed Lorentz gas in equilibrium the topo- 
logical pressure vanishes for (3=1. Now D((3) is given by 

D(f3) = 



d{v d + P (/3)) 



1 



(d-l)(l-(3) 



[l + (d-l)(l-(3)](3(3-d) 



(15) 



2[d- 1-1-/3(3 -d)\ 

From this we can easily get the normal diffusion coefficient 
by setting (3=1, i.e., 

v 2 (d+l) 



D(l) 



Adv d 



(16) 
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(a) 



(b) 
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FIG. 2: (a) Topological pressure divided by the collision frequency as 
a function of (3 for different dimensions d, with the parameter choice 
a = 1, v = 1, n = 0.001, and k = 0.001. (b) A close-up of this 
for /3-values around 1. 



FIG. 3: (a) Dynamical entropy divided by the collision frequency, 
as function of (3 for different dimensions d and with the parameter 
choice a = 1, v = 1, n = 0.001, and fco = 0.001. (b) A close-up of 
this for /3-values around 1. 



OTHER DYNAMICAL PROPERTIES 

As mentioned before, from the dynamical partition function 
we may obtain several other dynamical characteristics of our 
system. For (3 = 1 the topological pressure equals minus the 
escape rate 7, thus we have 7 = }CqD(1). From Eqs. (|4} and 
(|9} we obtain the dynamical entropy function h(j3) as 



h((3) = ho(J3) - k 2 [D{(3) - (3D'{(3)\ , 



(17) 



with D'{(3) = dD{(3)/d(3 and h ((3) the entropy function for 
the closed Lorentz gas in equilibrium. 

Figures 12 and |5] show the topological pressure and the dy- 
namical entropy, respectively, divided by the collision fre- 
quency, as functions of (3 for different dimensions d. As ex- 
pected, P((3) is negative for (3=1 because P(l) = —7, see 
the inset of Fig. [2 Furthermore, P((3) is a convex function for 
all dimensions considered. Since fco *C 1 for large systems, 
the deviations of P{(3) and h((3) from their equilibrium values 
are small. 

The logarithms of the prefactors of the correction terms in 
Eqs. and (II 71 proportional to k 2 are plotted in Fig.@] For 
P((3) this prefactor is D{(3), which is always positive within 
the allowed ranges of (3 and d. As can be seen from Eq. dl 51 
D{(3) diverges at a pole located at (3 = (d — l)/(d — 3). In 
Fig-IHb) the logarithm of the absolute value of the correction 
term for the dynamical entropy is plotted because the prefac- 
tor changes sign at (3 s» 0.4, see inset Fig. EJb). Thus the 
correction to the topological entropy at (3 = will be negative 
while the one for the KS-entropy at (3 = 1 will be positive 
(see also the discussion). 



The KS entropy for general d for the open Lorentz gas is 
given by 



Ml) = h KS = h° KS - k 2 D(l) 



1 



Vd 



(18) 



where h° K s is the equilibrium value of the KS entropy for the 
closed Lorentz gas in d dimensions. The specific form of this 
as a function of d reads 



n KS 



Vd{l ~ d) 



+^(3-d) 



C - In 



2v 
au d 
d 



(19) 



with C Euler's constant and fy(x) the digamma function 12011 . 

For d = 2 and d = 3 we can compare our results to pre- 
vious calculations based on an extended Lorentz-Boltzmann 
equation approach (kJ. From Eq. (1191 we get the KS entropy 

for d = 2 as 



l KS 



h° - h 2 



3_u_ 

8^ 



2 r 



h° 



V-2 



(20) 



h? KS is given by the one positive Lyapunov exponent in equi- 
librium 



h°KS — A — i/ 2 



1-C-lnf^ 
V v 



The KS entropy for d = 3 follows as 
hKS = h KS ~ k 07~- 



V3 



(21) 



(22) 
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FIG. 4: Natural logarithms of the correction terms proportional to kg 
of (a) the topological pressure, see Eq. and (b) the absolute value 
of the dynamical entropy, see Eq. <17> . as a function of (3 for different 
dimensions d. The inset in (b) shows a close-up of [D(/3) — (3D\(3)] 
around (3 = 0.4, where this function crosses zero. All results are for 
a = 1, v = 1, and n = 0.001 



Here, h° KS is given by the sum of the positive Lyapunov ex- 
ponents in equilibrium 







Kv) 


-c 







(23) 



As one should expect, Eqs. ( 1201231 are in agreement with the 
previous calculation. 

From the escape rate formalism l2lll we know that the KS 
entropy equals the sum of positive Lyapunov exponents minus 
the escape rate. Note that here the Lyapunov exponents are de- 
fined on the repeller. Furthermore, for (3=1 the topological 
pressure equals the escape rate 7 and this can be expressed for 
the open Lorentz gas as 7 = k^D{\). Therefore, we have an 
equation for the sum of positive Lyapunov exponents of the 
open Lorentz gas, 



i i 

~l-d 



klD(l) 



n KS 
Vd 



(24) 



(25) 



The correction term proportional to the diffusion coeffi- 
cient is always positive, because for low densities the term 




FIG. 5: Relative corrections to the topological and to the KS en- 
tropy, (a) shows the natural logarithm of (Titop — fttop) iy d/(7' l top) 
and (b) shows {has ~ ^Ks) v di '(jhjcg), both as functions of d, for 
parameter values a = 1, v = 1 and n = 0.001. 



h\(2v / avd) ^ 1 in Eq. d23i dominates. Hence, the sum of 
positive Lyapunov exponents is always larger for the open sys- 
tem than for the corresponding closed system. This again is in 
quantitative agreement with previous results Jlfill . 

Since we have an expression for general (3 we can also cal- 
culate the topological entropy h top which is given by h{(3) for 
(3 = 0. For general d this is given by 



Ho) 



Hop 



"top 



, 2 2D(l)u d 



h° 
'hop 



(26) 



with the equilibrium value of the topological entropy 



'hop 



Vd 



d-1 

~7T 



2v 
av d 



d-1 



(27) 

So we may conclude that the topological entropy is always 
smaller than in equilibrium. 

In Fig. |5] the relative corrections to the KS entropy and the 
topological entropy are plotted as functions of d. For h t0 p the 
ratio to the equilibrium value is very small. It decreases ex- 
ponentially with d, because asymptotically for large d h® op is 
independent of the radius a, whereas the collision frequency 
Vd is proportional to aSd — 1). For hxs tne correction, scaled 
as in the figure, appoaches unity in the limit of large d, as 
follows readily from eqs. Jl 8I > and © plus Stirling's approxi- 
mation In T(x) w x In x — x. 

More dynamical properties can be obtained from P{(3). 
The partial Hausdorff dimension is given by the value of (3 
where P((3) intersects the /3-axis, whereas the partial infor- 
mation dimension is given by the value of (3 where the tangent 
at P(l) intersects the /3-axis 14 12211 . For the latter we easily 
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find that the tangent is given by 

P(1)-(1-/3)P'(1) 

= (1 - (3)h° KS - k 2 [£>(!) 



(1-(3)D'(1)), (28) 



Thus, the partial information dimension follows as 



di = 1 - 



1 - 



k 2 D'(l) 
1 



(29) 



and is clearly smaller than unity. For a large system with fixed 
density an expansion for small escape rates 7 gives dj 1 — 
j/hjcg. The complexity of P(/3) and the /3-dependence of 
the diffusion coefficient prevent a calculation of the partial 
Hausdorff dimension. However, for large systems, where ko 
is very small, the partial Hausdorff dimension will be very 
close to dj because the point of intersection of P{(3) with the 
/3-axis will be close to (3 = 1, therefore in Eq. l|9} both terms 
are well approximated by a Taylor expansion around (3 — 1 
up to linear order in (3. 



DISCUSSION 

We may conclude first of all that for the Lorentz gas the cal- 
culation of corrections to the topological pressure due to open 
boundary conditions is remarkably simple; it just requires the 
solution of an effective /3-dependent diffusion equation with 
the same open boundary conditions. One question that could 
be asked here is whether the diffusion coefficient to be used 
in this calculation is the same as in a closed system. Since 
realizations of the random flight process with a slightly higher 
than average collision frequency will have a slightly lower dif- 
fusion coefficient these will lead to a slightly smaller escape 
rate from the system. The same can be argued for processes 
in which the frequency of backscattering events is slightly en- 
hanced and that of forward scattering is slightly suppressed 
with respect to the averages in an equilibrum system. There- 
fore the average diffusion coefficient on the repeller should 
be slightly smaller than the diffusion coefficent of the equi- 
librium system. However, one easily estimates that the sup- 
pression of the diffusion coefficient is of order k 2 , leading to 
a correction of order k^ in the escape rate. Our results for the 
corrections to entropies and dimensions, which are all of order 
k 2 , therefore will not be changed. For higher order corrections 
of course such terms are important. 

It is very interesting that for the dilute random Lorentz gas 
the effects of the open boundary can be separated straight- 
forwardly into effects due to changes in the distributions of 
free flight times and of scattering angles respectively. This is 
because the stretching factor first of all factorizes (to lead- 
ing order in the density) into a product of stretching fac- 
tors pertaining to a single free flight plus subsequent col- 
lision, and in addition each of those factorizes into factors 



depending on the free flight time respectively the scattering 
angles alone. Specifically, one has Ad(0,t) — A1A2, with 
Ai = (2vt/a) d ^ 1 and A 2 = cos d ~ 3 9. Rewriting the factor 
A 1- ' 3 as A] -131 A^~^ 2 one may obtain the topological pres- 
sure and the effective diffusion constant as functions of both 
(3i and (3%. By taking the derivatives of the topological pres- 
sure with respect to (3\ and (32 one finds that the two terms 
on the right hand side of Eq. J19I may be assigned to the dis- 
tributions of free times and of scattering angles respectively. 
Similarly the correction to the sum of the positive Lyapunov 
exponents, given by k 2 D'(l), can be separated into a term due 
to the change in the free time distribution, which is of the form 



A t EA+ = klD{l){d- 1) 



In- 



2v 
v d a 



-C- 



1 



(30) 



and a term due to the change in the distribution of scattering 
angles, reading 



k 2 



(3-d) 



D(l) 



1 



1 



1 



(31) 

One might be tempted to think that the changes in the distri- 
bution function are due primarily to particles near the open 
boundary, which will only survive if their free flights keep 
them inside the system. This however is completely false. At 
any given instant the fraction of particles inside a layer of a 
few mean free paths near the boundary may be estimated as 
proportional to 1/R 2 , with R an estimate for the diameter of 
the system; the volume fraction covered by the boundary layer 
is proportional to 1/R and the density near the open boundary 
is also of order 1/R compared to the average density. In ad- 
dition, in order for a particle trajectory to be on the repeller, it 
has to remain inside the system forever after. For trajectories 
near the open boundary at the given instant, the probability 
for this to happen is another order 1/R smaller than for tra- 
jectories at large [16]. Therefore trajectories getting near the 
boundary at any time do not contribute at all to the order /cq . 
Rather, the deviations from the equilibrium distributions are 
caused by the fact that free flights or scattering angles that 
move a particle away from an open boundary lead to higher 
survival probability than ones that bring it closer to it. If one 
wants to know how, locally, the distribution of free times or 
scattering angles is changed, he would have to take recourse 
to the methods laid out in Ref. 11611 . It is fortunate though, 
that the reduction of the dynamical partition function to an 
effective diffusion problem allows one to calculate the KS- 
entropy and the like, at least up to order k 2 , without having to 
go through the details of this formalism. 

An interesting question is in how far the present method 
of calculating corrections to the topological pressure from the 
solution of a diffusion equation can be generalized to systems 
of many moving particles. Unfortunately this does not seem 
to be possible in any straightforward way, even for dilute sys- 
tems with hard interactions. The reason for this is the semi- 
dispersing character of the dynamics of these systems, inter- 
preted as billiards in a high dimensional phase space. Due 
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to this property the dynamical partition function does not ap- 
proximately factorize into a product of terms describing the 
effects of one free flight and the subsequent collision. Note 
that even the Lorentz gas loses this property at higher densi- 
ties, where the mean free path between collisions is not large 
compared to the scatterer radius. 

Finally, we remark that, like in the equilibrium and field 
driven disordered Lorentz gas, the calculations of the topo- 
logical pressure for general values of j3 have to be taken with 
caution Points in phase space are weighted by 

A 1 "' 3 , thus for (3 sufficiently < 1 the dynamical partition 
function will be dominated by the largest stretching factors, 
which are due to trajectories confined to regions of high den- 
sities of scatterers. Even though the number of such trajecto- 
ries decreases exponentially with time, the stretching factors 
raised to the power 1 — p of the remaining ones will still make 
these dominant for p sufficiently far from 1. For /3 sufficiently 
> 1 on the other hand the dynamical partition function will be 
dominated by trajectories confined to the neighborhood of the 
least unstable periodic orbit. 

In conclusion, we have shown how to relate the thermo- 
dynamic formalism for the open Lorentz gas to a diffusive 
random flight problem. We have calculated the topological 
pressure in d dimensions as a function of the temperature-like 
parameter p. For the open Lorentz gas, the topological pres- 
sure is the sum of the topological pressure for the closed sys- 
tem and a /3-dependent effective escape rate which is given by 
a /^-dependent diffusion coefficient multiplied by the square 
of the smallest wave number fco fitting the diffusion equation 
with absorbing boundary equations. From this we have ob- 
tained several dynamical quantities such as the Kolmogorov- 
Sinai entropy, the topological entropy, and the partial infor- 
mation dimension for general dimension d. 

This work was supported by the Collective and coopera- 
tive statistical physics phenomena program of FOM (Funda- 
menteel Onderzoek der Materie). 
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